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Abstract. Nuclear Star Clusters are observed at the center of many galaxies. In 
particular, in the Milky Way's center, the Nuclear Star Cluster coexists with a central 
supermassive black hole. The origin of these clusters is still unknown; a possible forma- 
tion mechanism is the decay of massive globular clusters driven inward to the galactic 
center by dynamical friction and their subsequent merging. By investigating this sce- 
nario by means of sophisticated N-body simulations we found that this process could 
lead to a final product which actually shows many of the observed features of the Milky 
Way Nuclear Star Cluster. 



1. Introduction 

Nuclear Star Clusters (NSCs) have been observed in many early-type galaxies and late- 
type spirals. The Milky Way (MW) and a handful of other galaxies contain both a 
NSC and a massive or even super massive black hole (MBH). NSCs are among the 
densest star clusters observed, with effective radii of a few pc and central luminosities 
up to ~ 10 7 L Q . The formation mechanism of NSCs is still unknown. Two competing 
models are possible. In the gas model a NSC ca n form from the gas th at migrates to 
the center of the galaxy, where then forms stars (ISchinnerer et al.ll2008l) . Due to the 
inherent complexity of gas dynamics, the gas model is almost qualitative and difficult 
to study. Alternatively, in the merger model, massive clusters decay to the center via 
dynamical friction and merge to form a dense nucleus (Tremaine et al. 1975; Capuzzo- 
Dolcetta 1993b . Observations of NSCs in dE galaxies suggest that the majority, but 



not all, dE nuclei coul d be the result o f packing mass in form of orbitally decayed 
globular clusters (GCs) (ILotz et al. 2004). Numerical simulations have also shown that 
the basic properties of NSCs, including their shape, mass density profile, and mass- 
radius relat ion, are well reproduced in the merger model under a variet y of explored 
conditions (ICapuzzo-Dolcetta & Miocchill2008allbl ; iHartmann et alll201lb . 



2. The Galactic NSC: formation and evolution 



The MW NSC has an estimated mass of ~ 10 7 M Q (lLaunhardt et alj|2002h. and it hosts 
an MBH whose mass, ~ 4.3 x 10 6 M Q , is uniquely well determined iGillessen et al. 
d2009h . Since the relaxation time at SgrA*'s influence radius is robustly estimated to 
be 20 - 30 Gvr dMerrittll20Toh. the NSC has not had enough time to settle on the Bahcall 
& Wolf d 19761) (BW) cusp stable state. Indeed, the late-type stellar population exhibits 
not a cusp but, rather, a space density which i s flat or even falling toward SgrA*, inside 
a core of radius 0.5 pc (Buchhol z et al.l l2009). 
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Figure 1. Left panel: Spatial profile of the central NSC after 3, 6, 9 and 12 merg- 
ers. The dashed line is the fit to the NSC profile obtained at the end of the last 
merging event using a broken power law model. Right panel: Spatial profile of the 
central galactic region at regular intervals of time after the final merging. Time grows 
upward in the figure. 



2.1. The simulations 

For the first time, we performed self-consistent Af-body simulations to test the merger 
model for the formation of the MW NSC. These simulations, unlike in previous works, 
involve the presence of a n MBH at the ce nter of the Galaxy. The details of these sim- 



ulations are d escribed i n Antonini et al. | (20 11) and the most significant results are 
summarized in lCapuzzo-Dolcetta et all (1201 lab . In the simulations, the central MBH is 
represented as a point mass particle initially located at the origin of coordinates. The 
N-body galactic environment is sampled from a shallow density cusp (p oc r' 1 ^ 2 ), while 
the individual GC N-body representation has been done starting from a tidally limited, 
massive and compact King model. Thus, we performed full N-body simulations to 
study the consecutive infall and merging of a set of 12 GCs each starting from a galac- 
tocentric distance of 20 pc on an initially circular orbit. This small initial galactocentric 
distance is the result of previous GC orbital decay due to dynamical friction on a larger 
spatial scale. After the first GC decayed to the galactic centre, we waited for the NSC to 
reach a quasi-steady state and we added to the system a second GC. This procedure is it- 
erated until 12 clusters merged in the inner regions of the galaxy. Thus, the NSC is built 
up through consecutive infalls of GCs and, after its formation, we followed the long- 
term evoluti on of the NSC. Simulations are done with both th e high precision, p arallel, 



NBSymple (ICapuzzo-Dolcetta et al.ll2011bl) and PhiGRAPE (Harf st et al.ll2007h codes 



running on the CPU+NVIDIA TESLA C2050 platform in Roma, Sapienza (Italy), and 
on the Rochester Inst, of Technology (USA) GRAPE cluster. 



3. Results 

The merging occurs rather quickly: after about 20 crossing times the resulting system 
attains a quasi equilibrium configuration. The orbital energy loss leads to a full merge 
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Figure 2. Left panel: Break radius as a function of time (in units of the f rj i n fl) in 
the evolving model (empty circles). The solid line represents the best exponential fit 
to the data. The region limited by the two dot- dashed l ines c orresponds to the two 
estimates of the MW NSC core radius given in iMerrit j d2010t) . Right panel: Radial 
profiles of the anisotropy parameter f3 = 1 - (cr 2 , + cr 2 .) / (2cr 2 ) during relaxation 
displayed at regular interval of times. Thicker curves correspond to older ages. 



of the all set of decaying satellites which settle on a meta-stable state, whose following 
evolution is mainly due to collisional relaxation. The left panel of Figure Q] shows the 
spatial profile of the system after the complete merging of 3, 6, 9 and 12 clusters. As 
inferred from the best broken power law fit (dashed line in the left panel of Figure [T|), 
the merger product shows a density profile with a central nearly flat core and a slope 
out of the core reg ion, p <x r~ 18 , in good agreement with what is observed in the MW 
(Halle7etaD[l996). On the other hand, t he core radius o f the merger remnant is 4 times 
larger than the observed value (~ 0.5 pc. lBuchholz et al.ll2009h . Anyway we expect that, 
as a natural product of two-body relaxation, such initial parsec-scale core shrinks with 
time approaching the BW steady state in a time which is of th e order of th e relaxation 
time at the influence radius (r; n fl) of the central MBH, f r j n fl (IMenitti 120 1 00 . Thus, we 
followed numerically the collisional relaxation process of the system and found that 
it undergoes a self-similar evolution which keeps almost unaltered the external radial 
slope (~ -1.8) while the core shrinks reducing its size (see right panel of Figure [J). 
Actually, the evolution of the inner density plateau size of the NSC is shown by the left 
panel of Figure |2l which gives the break radius, 77,, of the best fitting broken power law 
profile as a function of time. The shrinking time dependence of the break (i.e. core) 

-2 24 i 

radius is well described by the exponential law r\,{t) - 1.57 e Right panel of 

Fig. |2] illustrates the NSC evolution toward isotropy by plotting the model anisotropy 
profile at different times. At the end of the simulation the NSC has only a small bias 
toward tangential motion. This is co nsistent with proper-motion data that indicate a 
slight degree of tangential anisotropy (Scho del et aLll2009klMerrittll20l0l) . 
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4. Conclusions 

Various authors inves tigated the merger model for the formation of Nuclear Star Clus- 
ters (see for example iTremaine et al.ll 19751 : ICapuzzo-Dolcetta & Miocchill2008alibh but 



in all these works the presence of a central MBH was neglected and the simulations 
were not finalized to a specific study of the MW. Our work covers this lack through 
self-consistent simulations whose ini tial conditions are set up basing on recent data 
about the MW dLaunhardt et al.ll2002h . including the presence of its central MBH. The 
merging of GCs produces a system which rapidly settles on a quasi-steady state, slowly 
evolving due to internal relaxation which, in its turn, is affected by the presence of the 
MBH. The central and almost flat core characterizing the density profile of the merger 
product is maintained for a time long enough to justify the core actually observed in 
the MW NSC. Our results are also supported by the fact that the external slope of the 
density profile of the NSC remains quite unaltered after the end of the last merger event 
and its value, ~ -1.8, is in optimal agreement with the one inferred from the obser- 
vations of the Galactic NSC. Moreover, we observe a drift toward velocity isotropy 
although keeping a slight degree of tangential anisotropy; this peculiarity was found 
in the observational data and already seen in the numerical studv bv Capuzzo Dolcetta 
et al. (I200lh . As a conclusion, the process of consecutive infalls and merging in the 
inner region of the MW of a number of the order of 10 massive globular clusters could 
give an explanation of the observed features of the NSC observed in our Galaxy around 
the Sgr A* massive black hole. 
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